TPI-MINN-13/97-T 
UMN-TH-1542/97 



An Illustrative Example of How Quark-Hadron Duality 

Might Work 

B. BLOK^ M. ShIFMAN, ^ AND Da-Xin Zhang^ 

^ Department of Physics, Israel Institute of Technology, Technion, Haifa 32000, 

Israel 

^ Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 54555, 

USA 



Abstract 

We discuss the issue of the local quark-hadron duality at high energies in 
two- and four-dimensional QCD. A mechanism of the dynamical realization of 
the quark-hadron duality in two-dimensional QCD in the limit of large number 
of colors, Nc — > oo, (the 't Hooft model) is considered. A similar mechanism 
of dynamical smearing may be relevant in four-dimensional QCD. Although 
particular details of our results are model-dependent, the general features of 
the duality implementation conjectured previously get further support. 



1 Introduction 



In the recent years the focus of apphcations of the operator product expansion (OPE) 
^ has shifted towards the processes with the essentially Minkowskean kinematics. 
Perhaps, the most well-known example is the theory of inclusive decays of heavy 
flavors (for a review see e.g. Ref. [§]). This fact, as well as the increased demand 
for more accurate predictions, puts forward the study of the quark-hadron duality 
as an urgent task. 

A detailed definition of the procedure which goes under the name of the quark- 
hadron duality (a key element of every calculation referring to Minkowskean quan- 
tities) was given in Refs. @, In a nut shell, a truncated OPE is analytically 
continued, term by term, from the Euclidean to the Minkowski domain 0. A smooth 
quark curve obtained in this way is supposed to coincide at high energies (energy 
releases) with the actual hadronic cross section. 

If duality is formulated in this way, it is perfectly obvious that at finite en- 
ergies deviations from duality must exist. The difference between the measured 
physical cross section and a smooth OPE prediction will be referred to as an os- 
cillating/exponential (duality violating) component. In Ref. it was shown that 
if we knew the leading asymptotic behavior of the high order terms in the power 
series we could evaluate this component. Unfortunately, very little is known about 
this aspect of OPE, and we have to approach the problem from the other side - 
either by modeling the phenomenon |^ or by studying some general features of the 
appropriate spectral densities. One can also try to approach the problem purely 
phenomenologically (for recent attempts see e.g. P). 

An illustrative spectral density, quite instructive in the studies of the issue of 
the quark-hadron duality, was suggested in 0, 

TV °° 

Im n = Const. ^) (1) 

2 n=l 

where 

and from now on we will drop an inessential constant in front of the sum. The color 
factor Nc is singled out for convenience. The imaginary part above represents, for 
positive values of s, a sum of infinitely narrow equidistant resonances, with equal 
residues. The distance between the resonances is A^. It defines n(g^) everywhere 
in the complex plane g^, through the standard dispersion relation, up to an addi- 
tive constant which can be adjusted arbitrarily. It is not difficult to see that the 
corresponding correlation function 



Ztc 



(2) 



^Moreover, usually one deals with the practical version Q of OPE, see Ref. js) for further 
details. 
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where ip is the logarithmic derivative of Euler's gamma function, and 




-S. 



In the Minkowski domain 8 is positive, in the Euchdean domain e is positive. Then, 
the asymptotic expansion of n(g^) in deep Euchdean domain is well known. 



where Bn are the Bernoulli numbers. At large n they grow factorially, as Bn ~ {2n) ! 
(see 0, page 23). We deal with the sign alternating series. 

Although the spectral density (|l|) is admittedly a model, it was argued that 
a similar factorial growth of the coefficients in the power (condensate) series is a 
general feature. 

The spectral density (|l]) may be relevant in the limit of the large number of 
colors, Nc —>■ oo, when all mesons are infinitely narrow. This limit is not realistic, 
however. Moreover, in this limit the local quark-hadron duality, as we defined it, 
never takes place since even at high energies the hadronic spectral density never 
becomes smooth, even approximately. One can smear it by hand, of course, but 
then deviations from the local duality will be determined not only by the intrinsic 
hadronic dynamics, as is the case in the real world, but also by particular smearing 
procedure - they will depend on the weight function chosen for smearing, the interval 
of smearing and so on. In the actual world the smearing occurs dynamically, since 
at high energies the resonance widths become non-negligible. The limits of — > cxd 
and iVc — >■ oo are not interchangeable. 

Here we suggest and study more realistic (dynamically smeared) spectral densi- 
ties compatible with all general properties of Quantum Chromodynamics. Starting 
from infinitely narrow resonances, as in Eq. (|l]), we introduce finite widths, ensuring 
smooth behavior. Technically, in the first part of the paper the problem of duality 
is analyzed in the two-dimensional 't Hooft model |^ (see also - [jl^]). The quark 
confinement in this model is built-in. We then try to abstract general features of 
this solution, which may persist in QCD. In the second part of the paper an attempt 
is made to work out the same mechanism in four dimensions. 

Dynamical "smearing" of the spectral densities taking part in QCD due to non- 
vanishing resonance widths, exhibits the same features of the high-energy behavior 
as was suggested in Refs. [2] and [3] on the basis of rather naive models, e.g. instan- 
ton models. Namely, the approach of the spectral density to the smooth limit (the 
deviation from duality) is exponential, with oscillations. This pattern seems to be 
general and may be considered now as a well-established model-independent fact. 
At the same time, the exponent determining the rate of the damping of the duality 
violations depends on details of the large distance dynamics. We were unable to 
find it from first principles, and had to settle for model-dependent determinations. 



n(g^)--§ ine + ^-E(-i) 



-Bn£_ 

2n 



-2n' 



n-1 



(3) 
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Let us note that several useful results on the relation between duality violations 
and the divergence of the power expansion were presented in Ref. [|^, and we 



incorporate them H. The very idea of using the 't Hooft model as a theoretical 
laboratory adequate to the problem was formulated there. Moreover, it was noted 
| |r3| that the leading asymptotics of the high-order coefficients does not require the 
knowledge of the exact mass spectrum. Suffice it to know the leading term in the n 
dependence, where n is the radial excitation number, which immediately translates 
in the leading factorial behavior of the coefficients. In particular, the position of the 
low-lying states is irrelevant in the regime with the factorially growing coefficients. 
This is a specific feature of the asymptotic series. 

The if) function model was originally suggested in Ref. [0 for the heavy-light 
quark systems. In Ref. |T3] it was noted that it was more appropriate for the light- 
light quark systems since in the heavy-light systems the resonances are not expected 
to be equidistant. Straightforward quasi classical estimates yield in this case that 
the meson energies (measured from the heavy quark mass) asymptotically scale as 
\/n. (See also Ref. [Q, where this scaling law is reproduced in a linear potential 
model.) In the present paper we further develop the if) function ansatz adapting it 
for the light-quark systems. 

The paper is organized as follows. Section 2 briefly reminds the basic strategy 
of the OPE-based calculations. Here we outline the main elements of the analysis 
to be opresented below. In Sect. 3 we flrst discuss the 't Hooft model in the 
leading 1/Nc approximation. Then we calculate the resonance widths and evaluate 
the polarization operator of the two scalar currents in taking into account 1/Nc 
corrections. In Sect. 4 we present qualitative arguments why the implementation of 
duality obtained in the 't Hooft model must be also relevant, at a qualitative level, 
in four-dimensional QCD. Sect. 5 gives conclusions and outlines possible physical 
applications of our results. 



2 OPE-Based Strategy, Sum Rules and Duality 

The issue of "deviations from duality" caused much confusion in the recent literature. 
Therefore, to begin with, we remind the basic strategy revealing, in clear terms, 
where the deviations can occur and where there can be no deviations. 

First of all, all calculations based on the operator product expansion are carried 
out in the Euclidean domain. Only away from the physical cuts this procedure 
is well-deflned. Any calculation consists of several crucial elements: identiflcation 
of the operators which can appear in the expansion, separation of hard virtual 
momenta (higher than a normalization point /i) from soft virtual momenta (lower 
than /i) and, flnally, calculation of the expansion coefficients in front of relevant 
operators. The latter is carried out in terms of quarks and gluons. That is why 

^We strongly disagree, however, with some particular calculations and expressions, presented 
in this paper. 
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the normalization point /i must be chosen sufficiently high, and the calculation 
must be done in the Euclidean domain. Even if in some problems calculations are 
conveniently presented in such a way as if they were done in the Minkowski domain, 
actually the corresponding results must be understood as an analytic continuation. 

The connection between the Euclidean predictions and measurable quantities 
is established via dispersion relations. In this way one can get certain sum rules. 
A large variety of them is offered on the market. If the Euclidean quantity is 
appropriately chosen, the (Euclidean) operator product expansion converges. The 
best-known example of such an appropriate choice is provided by the SVZ sum rules 
T5| obtained by virtue of the Borel transformation of the dispersion representation. 



Consider, for instance, the model of Ref. 0, see Eq. (0). The l/Q^ expansion 
of n(g^) in the Euclidean domain is factorially divergent. At the same time, the 
expansion of the Borel-transformed quantity has a finite radius of convergence. 
Indeed, the SVZ sum rule in the case at hand has the form 



iV,A2 e-'^VM^ 1 r 



(4) 



(in the right-hand side A^ plays the role of the inverse slope of the Regge trajec- 
tory). Assume that the left-hand side was calculated theoretically, using OPE, as 
an expansion in The domain of convergence of the power expansion is de- 

termined by the position of the nearest singularity in the complex plane. It is 
quite obvious that the expansion converges at 

ZTT 

Not only the radius of convergence is finite, due to the factor 2tt in the denominator 
the domain of convergence extends to quite low values of M^. If A^ ~ 1 to 2 GeV^, 
the power series is convergent at as low values of as 0.4 GeV^. This fact was 



empirically observed long ago |15 



If we consider the SVZ sum rule inside the convergence domain of OPE (see 
Eq. (H) , and are aimed at predicting the exponential integral on the right-hand side 
of Eq. (^) per se, it is meaningless to speak about deviations from OPE. In this 
formulation of the problem there are no deviations. 

The problem of deviations arises when we try to predict the spectral density Im 
n point-by point, at large s, or certain integrals of Im 11, not directly reducible to 
"good" sum rules of the type (^). Certainly, if we assume that Im 11 is smooth start- 
ing from some boundary value sq, Eq. (^) allows us to predict Im n(s) at s > sq 
unambiguously. If an oscillating component is allowed, however, one can always 
invent such a wild oscillating function, which, being integrated with the exponential 
weight, gives a contribution on the left-hand side of the sum rule (^ less than the 
last term of the power expansion retained, no matter how small this last term is. 
This component is referred to as duality-violating. Clearly, Quantum Chromody- 
namics admits only very specific oscillating/exponential components in the spectral 
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densities, if at all. The question is what particular oscillating/exponential functions 
are allowed by QCD dynamics. We first try to answer this question in a simpler 
context of the 't Hooft model, and then pass to discussion of QCD. 



3 The 't Hooft model 

The aim of this section is to study the quark-hadron duality in the 't Hooft model. 
We shall first briefly review the 't Hooft model in the Nc ^ oo limit and discuss 
1/Nc corrections in this model. Our original contribution is calculating the resonance 
widths for high excitations. Then, using these widths, we suggest an ansatz for the 
asymptotic behavior of the polarization operator of two scalar currents, based on 
the Breit-Wigner approximation. This resonance-saturated polarization operator 
will be referred to as phenomenological. We will confront it with the truncated 
power expansion. The difference between these two expressions gives an idea of the 
duality violation. 



3.1 The 't Hooft model: generalities 

Two-dimensional QCD is described by the Lagrangian 

L = -\G%G% + q^iiP-mf)q^ , (6) 

where / is the flavor index. Since the multifiavor aspect is irrelevant for our problem, 
we shall consider, for simplicity, one flavor; correspondingly, the index / will be 
omitted hereafter. If the gauge coupling of the theory is it is convenient to 
introduce an effective coupling, 

f = g'N,. (7) 

which stays constant in the limit N^. oo. The coupling constant g has dimension 
of mass; it sets the scale for all dimensional quantities in the chiral limit, i.e. m -h> 0. 
We introduce the scale 

^ = — = —r ' ^) 

and measure all quantities in these units, e.g. the quark mass 

2 2/2 

while the mass of the n-th meson 

2 2/2 

and so on H. 



■^Normalization of the coupling constant g coincides with that adopted in Refs. and jo), but 
differs from that in Ref. 
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As was shown by 't Hooft, the model with the Lagrangian (^) is exactly solvable 
in the limit Nc ^ oo. The bound state spectrum includes an infinite number of 
bound states whose masses lie on an almost linear trajectory. The properties of 
these bound states are described by the 't Hooft equation, 

_ (7' - i)0n(x) My)dy 

'^"^"^ ^" x(l-x) Jo {x-yf ^ ^ 

Here the integral is understood as a "principal value" , 



— = lim - 

p2 e^O 2 



1 1 

+ 



{p + ieY {p — it) 



Moreover, x is the momentum fraction of the meson carried by the antiquark (in 
the infinite momentum frame), while 1 — x is that of the quark; 0„(a;) is the wave 
function of the n-th bound state. 

The integral equation (^ must be solved with the boundary conditions: 

- { (1 1: 1 1 . (10) 

Here /3 is the smallest (in the absolute value) root of the equation 

7r/?cot(7r/?) = 1 - 7^. (11) 

Below we shall be interested in the massless case, 7 = 0. In the massless limit, 
the lowest-lying state (the "pion") can be found from Eq. (|^) analytically. Indeed, 
the solution 0o = Const., corresponding to /io = 0, obviously goes through. We are 
interested, however, in highly excited states, 1. 

In the original paper 't Hooft suggested the following approximation for 0„(x) 
at large n, x not too close to and 1: 



[X 



sin(n7rx). (12) 



Recent calculations |]T6|, exploiting a new improved numerical procedure (the so 
called spline method), show that for the massless case a better approximation for 
large n is 

(j)n{x) = V2cos{TTnx) . (13) 

This formula works very well numerically everywhere except the very endpoints 
X = 0,1, and a slight x-dependent shift when x ~ 0.5. 

Note that the wave functions (pTsD satisfy the proper boundary conditions for the 
massless case, which corresponds to /? equal to zero, and, hence, 

0„(O) = C,0„(1) = PC. (14) 

Here P = 1 for the states with the even parity and P = —1 for the states with the 
odd parity, and the constant C 7^ 0. 
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The mass spectrum in the massless case was found by 't Hooft. The asymptotic 
behavior of /i,^ is 

= TT^n [1 + 0(ln(n)/n + ...] . (15) 

Note that the mass formula (|15|) does not depend on the choice of the wave 
function, Eq. (|13D or the 't Hooft choice, at least in the leading in n approximation 



3.2 The meson widths 

As was already mentioned, in the limit Nc —>■ oo the bound states in the 't Hooft 
model are stable, their widths vanish. However, once one takes into account the 
leading 1/Nc correction, the resonances begin to decay. In the first order in 1/Nc 
expansion there are only two-particle decays a —>■ b + c. The relevant coupling 



constants gate are given by the following formula P, O, 16, 10, 11 



9a,c = f^'J^il - (-l)('^"^'^^+'^^)](/l + faJ ■ (16) 



c 



Here aa is the parity of the a-th resonance. The constants are determined from 
the following expressions: 

2 r^ziz fx uj-i- \ 

fabc = / Mx)M^/^±)^c 



UJ± Jo \^ — ^± J 

-—(' 0.(x)<l>,(xM)0, f ^^"l , (17) 

UJ± Joj± \1 ~ UJ± J 

where u± are two roots of the algebraic equation corresponding to the mass-shell 
condition, 

ml = ^ + j^y (18) 

00 [l — iO) 



The function ^a{x) is defined as 

Jo [x — 



My) 



Using the above expressions for the decay couplings one can readily calculate the 
resonance widths in the leading 1/Nc approximation. They are given by 

where / is the standard "triangular" function, 

1 

/(m„, nib, rric) = -[ml - {nib + m^f] [ml - ijUb - mj^] . (20) 



The sum in Eq. ( |T9| ) runs over all mesons b and c with the constraint mb + mc < ma- 
Our task was to establish the asymptotic behavior of the widths, as a function of 
the excitation number, at large values of n, in the leading 1/Nc approximation (all 
widths are proportional to 1/Nc] the excitation number n will be temporarily called 
a in this section). In order to find the widths we first computed analytically, using 
the wave functions (p!3D, the overlap integrals (|1^). The answer can be expressed via 
the integral sine and cosine functions and was obtained using the REDUCE program. 
Since it is very bulky it seems unreasonable to present here the final expression 



After computing the overlap integral we performed numerically summation over all 
possible b and c in Eq. ( pj] ) for a up to 500. The result for the widths exhibits 
a remarkable pattern. The widths of the individual levels oscillate near a smooth 
square- root curve, see Fig. 1. This figure shows the width of the a-th state vs. a, 
up to a = 500. The result of averaging over the interval of 20 resonances is depicted 
in Fig. 2. We see that the curve of the averaged resonance widths r(a) = Ta is very 
well approximated by the function 

T{a) = 4^V^{l + 0{l/a)), (21) 

where the parameter n is introduced in Eq. (H), and A is a constant which will be 
given below. 



Since the square-root law (21) for the (averaged) widths is valid in such a large 
interval of the excitation numbers and turns out to be so accurate, it seems plausible 
that this formula could be obtained analytically. This is an interesting question by 
itself, especially in four-dimensional QCD. Unfortunately, we were unable to find 
exact analytic solution so far. Some qualitative arguments in favor of the exact 
square root dependence are discussed in Sect. 4. The numerical value of the constant 
A is 

A = 0.44 ±0.05. (22) 

Below this result for the (averaged) widths will be used for determining the asymp- 
totic behavior of the polarization operator 0. 

Since the 1/Nc result for the decay width grows with a, one may worry about the 
1/iV^ corrections. If they grew with a sufficiently fast this could invalidate Eq. (pT]) 
in the interval of the excitation numbers we are interested in, namely, a = Const. 
Nc, where the constant above can be numerically large, but it does not scale with 
Nc- Using quasiclassical arguments (see Sect. 4) one can show that the actual l/N, 
expansion parameter in Eq. (^) is ^/a/Nc. This means that at a = Const. Nc 
corrections 0{1/N^) and higher are negligible. 

Concluding this section let us note that the same square-root was reported pre- 



c 



viously in Ref. [0. We failed to reproduce the arguments of this work leading 



^Let us note that if one calculates the widths with the wave functions (p^, the square-root 
behavior of Eq. ( pl| ) is intact, but the value of the constant A is different, A' ~ 0.007, i.e. ~ 50 
times smaller. This fact indicates that the square-root law is not sensitive to the precise form of 
the wave functions, while the value of the coefficient A is. 
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Figure 1: The width of the a-th state (in units of /in ^N^ ^) as a function of a up 
to a = 500. 
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Figure 2: Smeared widths in the same interval of a and the square-root fit, Eq. (21). 
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to the square-root law, however. What is important, the constant analogous to 
A in Ref. fl^ is claimed to be proportional to l/y/m (!), and, thus, blows up for 
massless quarks. This poses perplexing questions. The coincidence looks completely 
accidental. 



3.3 The Breit-Wigner approximation. An Ansatz for the 
polarization operator 

Once we had found the resonance widths, we can calculate the polarization operator 
in the Breit-Wigner approximation. Here we shall consider the most interesting 
case of the polarization operator of the two scalar currents. Let us start from this 
polarization operator in the Nc —>■ oo limit |^ |19 . 



n(?^) = - E , . (24) 



Define the two-point function of the scalar currents j = qq, 

U{q') = t j rfV^(0|T{j(a:), j(0)}|0) . (23) 
In the Ai"c ^ oo limit 11 is given by |p 

i 

Here the constants Qn are the current residues 

{0\j{x)\n)=g^. (25) 

Note that the residues vanish for even n [0, so that the sum in Eq. (|2^) runs over odd 
n only. Below we shall be interested in the behavior of H(g^) for large ^ /i^. This 
behavior is dominated by the terms in the sum (plj) with large n ^ 1 |^, [13| . In 
order to calculate this sum explicitly we then need the large n behavior of Qn- It was 
determined in the same classical paper [H from the requirement of the compatibility 



of the expansion ( p^ and the perturbation theory asymptotics in the — > oo limit, 

H(g^)^-^ln^, Q'^-q\ (26) 



The coefficients Qn for sufficiently large odd n must be independent of n and are 
equal to 

9n = NcTTfx'. (27) 

Then, taking into account the linear n dependence of mass squared one can approx- 
imate the polarization operator (^3]) for sufficiently large |g^| by the ip function. 



n(g2)-H(0) = -^^(cr), a = ^ + -. (28) 



We hasten to emphasize again that this formula is not supposed to work at non 
asymptotic values of Q^. For instance, it does not contain the massless "pion" 
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Moreover, by shifting a little bit the masses and residues of the low-lying resonances 
we let n(g^) "breathe" at small without changing the asymptotic behavior. 

What will happen if we take into account finite widths of the resonances? To 
answer this question we calculate n(Q^) in the Breit-Wigner approximation. (Con- 
tinuation of the Breit-Wigner formula in the complex plane, away from the resonance 
position, is not unambiguous. We choose a specific continuation leading to proper 
analytic properties of the polarization operator, see below.) The inverse propagator 
of the ra-th bound state can be written as 

D-\q') = -{<f~< + nq')) (29) 

where S(g^) a function of order 1/Nc reflecting the possibility of the transitions 
a ^ be ^ a. This function is known at = m^, 

/I (TOO 2 

Im S(g2 = ml) = m„r„ = — -f . (30) 
Here we used Eqs. (13) and (19). Now we can write 

^n' = Q'(l-^l^|j)+"^n- (31) 

It is easy to see that (i) at = Eq. (|30|) is satisfied; (ii) the pole is shifted to 
an unphysical sheet, so that on the physical sheet there are no singularities except 
the cut at positive real q^. The property (i) is quite obvious. Let us comment on 
the property (ii). 

The easiest way to demonstrate that there are no singularities on the physical 
sheet is as follows. Observe, that at our level of accuracy one can write, instead of 
Eq. dTD, 

D-' = {z + ml) (32) 

where 

z = Q2(gV27rV)'^ • (33) 

Here the constant in Eq. (^) is adjusted in accordance with Eq. ( |35[ ) below. 

The physical sheet on the complex plane (Fig. 3a) is mapped onto a sheet 
with a "defect angle" on the complex z plane (Fig. 36). Going into the shaded area 
we pass to the unphysical sheets. Note that the pole of Dn lies in the shaded area. 

Assembling all pieces together we conclude that, with the resonance widths 
switched on, Eqs. (p3) and (E3) are substituted by 



H(g2) - n(0) = Const. Y: D-\Q') = - ^ _ Ajin^N^) £ ^^^^ 

where 

a = + - . (35) 
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Lower side of the cut aIhi z 



ReQ 




Rez 



Upper side of tlie cut 



(a) 



(b) 



Figure 3: Analytical structure of the polarization operator, (a) The polarization 
operator must be analytic everywhere in the complex plane, except the cut 
running on the negative real semi-axis of (positive real semi-axis of g^). The 
imaginary part of the polarization operator must be positive at the upper side of 
the cut; (6) The mapping of the plane onto the z plane, Eq. (33). The physical 
sheet on the Q"^ plane corresponds to the z plane with the shaded sector removed. 
The boundaries of the sector correspond to the lower and upper sides of the cut. 
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The constant in front of ipia) is adjusted in such a way as to leave intact the high 
asymptotics, see Eq. (^). It is clearly subleading in l/Nc- Strictly speaking, we 

should have omitted it at the level of accuracy accepted here. 

By construction, all singularities of the polarization operator ( plf ) are on the 

unphysical sheet. The discontinuity at the cut > will be calculated below. 



3.4 OPE and the asymptotics of the polarization operator 

The expression (32) for the polarization operator is clearly not exact, since we have 
made a number of approximations. They do not affect, however, the large 
asymptotics of n((5^). Neither the leading asymptotics of high orders in the 
expansion depends on these approximations. Since we are interested in the high 
energy behavior, distortions introduced in Il{Q^) at finite by the approximations 
made are not important. In particular, we will disregard the fact that low-order 
terms in the 1/Q^ expansion of Il{Q^) ("condensates") may come out with "wrong" 
coefficients. By adjusting the positions and residues of a few lowest resonances we 
can always get any desirable coefficients for any given finite number of terms in the 
l/Q^ expansion. 

Let us first study the impact of the resonance widths on Il{Q^) in the Euclidean 
domain. As we saw, the effect due to the nonvanishing widths essentially reduces to 
the substitution of the variable a defined in Eq. (26) by a, see Eq. (33). Therefore, 
the change in the asymptotic l/Q^ expansion at large (Euclidean) values of is 
rather insignificant. If at iVc = oo the n-th term of the power expansion is CniQ"^)'^" , 
at the 1/Nc level it becomes 

^» (g2)Mi-^) ~" ^"(Q^^^ + afcnlnQ') , (36) 



« = ^- (37) 



A 

The second term on the right-hand side in Eq. ( P^ ) is a small correction which 
reminds a logarithmic anomalous dimension in QCD. Unlike the first term in Eq. 
(PB|), the logarithmic term develops an imaginary part at large positive real values 
of q^. If we treated Eq. (|36D in the framework of OPE (more exactly, in practical 
version |Q), then we would predict that the spectral density 



Imnfsl 



2 



2 



" = "0 fy^t. 



n=l \ ^) 



(3J 



where s = and no is the highest term retained in practical OPE. By construction, 
the prediction for the spectral density obtained from practical OPE is smooth. As 
a matter of fact, all correction terms in Eq. (|38D are suppressed by l/N^, since 
a ~ l/Ai'c, and are numerically insignificant at large A^^^. 
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Let us examine now the "physical spectral density", i.e. the imaginary part at 
positive following directly from Eq. (|3^). Our task is to reveal an oscillating 
component, not suppressed by l/iV^. 

In order to find the imaginary part analytically, it is convenient to use the re- 
flection property of the if) function, 

ip{a) = ip{—a) — Ticotijia) — Ija . (39) 



It is obvious that the polarization operator ([3^ ) is a sum of three terms, correspond- 
ing to three different terms on the right-hand side of Eq. (pUj). Moreover, it is 
evident that the first and the third terms are smooth functions of q^. Their contri- 
bution to the imaginary part of the polarization operator for large s corresponds to 
the smooth component obtained from OPE, see the correction terms in Eq. (|38|). 

In order to study the oscillating/exponential component one must consider the 
second term, 

T m ^ ^-T ,( ~\ sinh(27/) 

Im n s = — Im cot TTfT = — — — - , 40 

2 2 cosh(2?/) — cos(2x) 



where 



a; = TrRecT ~ — ^_ ^ (1 — Q;ln(s)) 



1/ = TTlmcr ^ -— . (41) 



27r^2 

as 

Taking into account only the leading in terms it is easy to rewrite the spectral 
density (^0|) as follows 



Imn(g^) ^ ^[1 + 2exp (-— )cos(— )] . (43) 



N sinh(^) 
'°-"M-f co.,(g)-;o,^) - («' 

At as <^ /i^ we are in the resonance zone (here and below we use the nomenclature of 
Ref. [3], see Sect. 5.2). In order to find the high energy behavior of the imaginary 
part of the polarization operator, we must go to the oscillation zone, to energies 
as ^ li'^. Then 

y[l + 2exp(--)cos(- 

The unit term corresponds to the leading asymptotics (it reproduces the OPE pre- 
diction), while the second term is an oscillating/exponential component, missing in 
practical OPE i. It presents a deviation from duality we are hunting for (remember, 
by duality we understand a specific procedure formulated in Refs. [|], H). 

The oscillating component is suppressed by the damping exponential, exp {—as/ /i^) 
A comment is in order regarding the exponent, 

as _ A s 



(44) 



^In terms of an index a introduced in Sect. 5.2 of Ref. Eq. ( ^3|) implies that ct = 2. 
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Figure 4: The spectral density corresponding to Eq. (34) in the 't Hooft model, 
versus s (in the units 27r^/i^). Nc is put equal to 3, and the normalization factor 
is chosen in such a way that asymptotically the spectral density dsiplayed must 
approach unity. 



This exponent determines the boundary between the oscillation and the resonance 
zones, So ~ n'^Ncfi'^. This estimate is in accord with intuition. Indeed, the exponen- 
tial suppression of the oscillations should start when the resonance width becomes 
larger than the distance between the neighboring resonances. Since r„ ~ and 
the distance between the neighboring resonances Am„ ~ nfi/y/n, this occurs at 
the excitations number n ~ n'^Nc. Thus, the factor 1/Nc in Eq. ( ^41) is obvious. 
What is more remarkable, is a large additional numerical suppression, vr^, pushing 
the boundary to higher energies and making the exponential suppression weaker. It 
can be traced back to the numerical suppression of the width in Eq. (|2T|) which, 
in turn, is due to a strong numerical suppression of the phase space. Because of 
this fact the damping of the oscillations occurs very slowly in the 't Hooft model. 
Figure 4 shows Im 11 in the window from the 35-th to 65-th oscillation. We see that 
even here the oscillation amplitude is quite sizable. Unlike other features, expected 
to be generic, this numerical aspect - a very slow rate of the oscillation damping - 
is specific to the 't Hooft model and is not expected to survive in four- dimensional 
QCD, since it is entirely due to "abnormally" narrow widths of the resonances in 
two-dimensional QCD. Semiquantitative estimates of the damping exponent in four- 
dimensional QCD will be given in Sect. 4. 

In summary, for < the polarization operator is given by the function of 
the positive real argument and is a smooth function with no singularities. It is well 
approximated by its asymptotic expansion. The latter is most conveniently writ- 
ten in terms of the variable a defined above (see Eq. (|35|)), for which it becomes 
the standard asymptotic expansion of the ip function. In the Minkowskean domain, 
starting from the scale sq ~ vr^yU^A^c, one can approximate the polarization oper- 
ator by its smooth asymptotics (analytically continued from the expansion in the 
Euclidean domain) plus an exponentially decreasing and oscillating term (^3]). 
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4 Four-Dimensional QCD 



Unlike the 't Hooft model, we cannot solve four-dimensional QCD. At the qualitative 
level, however, we can apply the same ideas. The general pattern of the asymptotic 
behavior will be the same, but the numerical aspects look different. 

Our goal in this section is to discuss the polarization operator of two vector 
currents: 

n^.(g') =^ / exp{tqx){0\T{j^ix)UO)}\0)d''x. (45) 
Here j^(x) is the vector current, 

jf,{x) = u-ff,d{x) . (46) 

Due to the conservation of the vector current the polarization operator can be rep- 
resented as 

n,.(g^) = (gV-gV)n(g'). (47) 

Furthermore, 11(5'^) is dimensionless, and is perfectly analogous to the polarization 
operator (21) we dealt with in the 't Hooft model. According to the standard 
wisdom of multicolor QCD, we expect that ai Nc = oo and high energies H(g^) is 
representable as a sum of equidistant infinitely narrow resonances, with a constant 
residue, i.e. we arrive at the same ip function. The dynamical smearing is again 
provided by the resonance widths (a 1/A^c effect). What is known about the widths 
of the highly excited states in four-dimensional QCD? 

In the case at hand we do not have in our disposal quantitative tools which would 
allow us to calculate the widths. Such a calculation could have been performed in a 
non-critical string theory, were this theory available. Unfortunately, the issue is not 
worked out i, and we have to resort to qualitative arguments. If the string-based 
picture of color confinement is indeed valid, one can hardly avoid the conclusion 
that the resonance widths must grow linearly with m„ [^, much in the same way 
as in the 't Hooft model. 



Let us remind quasiclassical arguments of Ref. showing that r„ ~ rUn/N, 



in four-dimensional QCD. When a highly excited meson state is created by a local 
source, it can be considered, quasiclassically, as a pair of (almost free) ultrarelativis- 
tic quarks; each of them with energy m„/2. These quarks are created at the origin, 
and then fly back-to-back, creating behind them a flux tube of the chromoelectric 
field. The length of the tube L ~ m„/A^ where A^ is the string tension. The decay 
probability is determined, to order 1/Nc, by the probability of creating an extra 
quark- antiquark pair. Since the pair creation can happen anywhere inside the flux 
tube, it is natural to expect that 

1 B 

Tn LA^ = — m„ (48) 



^Some attempts in this direction were reported in Ref. [pO[. 
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where 5 is a dimensionless coefficient of order one. 

Let us note in passing that the corrections due to creation of two quark 

pairs are of order L'^ /N^ within this picture. Since L ~ m„ ~ ^/n^ the expansion 
parameter is y/n/Nc. 

Strictly speaking, the estimate (HSf) must be, rather, viewed as a lower bound. 



since the transverse fluctuations of the flux tube can increase the decay probability. 
However, most likely, these transverse fluctuations will materialize in the form of 
emission of glueballs, a sub leading l/N"^ effect, which is not considered here. It 
is not fully clear what impact these fluctuations may have on the quark- ant iquark 
pair production. It seems plausible that they only affect the numerical coefficient 
in Eq. (^8|), which is not calculated anyway. Given all naivete of the arguments 
and the estimate P8|), at the present stage it is reasonable to accept it as a working 
hypothesis. 

As was mentioned, in real QCD we expect i? ~ 1, see estimates in Ref. and 



below. Apart from this numerical difference, everything else is perfectly analogous 
to the consideration we have carried out in the 't Hooft model. Hence, we conclude 
that in four-dimensional QCD a viable model of the approach to asymptotics is 
provided by the same 4) function 

H(Q2) _ n(0) = Const. %lj{a) , (49) 

where in the case at hand 

Here C is a numerical constant correlated with the position of the lowest resonance, 

ml = C^m' (1 + 0(l/iVe)) , ^ = I- (1 + 0{l/N,)) . 

mo Nc 

The normalization constant in front of the ip function in Eq. ( ^9]) will be chosen in 
the form 

-11-^ 



so that the imaginary part at asymptotically high s is automatically normalized to 
unity. The plot of Im H is shown on Fig. 5. Although fine details of this spectral 
density are somewhat distorted compared to the actual experimentally measured 
spectral densities in the e~^e~ annihilation or r decays (e.g. the p meson width 
comes out ~ 1.5 larger than the experimental one) the qualitative similarity of our 
model with experiment is remarkable. To substantiate the point we display a plot 
of the spectral density in the vector isoscalar channel borrowed from Ref. p2| (Fig. 
6). 
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Figure 5: The imaginary part of the polarization operator corresponding to Eq. (50) 
with Am^ = 1 GeV^, C = 0.6 and B = 0.78. The number of colors is set equal to 
three. The horizontal axis presents s, in GeV^. 



P 




0.5 1 1.5 2 2.5 3 3.5 

Mass^ (GeW/c^f 



Figure 6: The spectral density in the vector isoscalar channel measured in e^e 
annihilation and r hadronic decays. 
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The oscillation zone starts at Sq ~ Am'^Nc{2TTB) ^. Above this boundary the 
asymptotic form of the oscillating exponential terms, analogous to Eq. (|^), is given 
by 

Imn(Q'^) = Const, [l + 2exp (-^) cos {^)] . (51) 

Both, the onset of the oscillation zone and the oscillation structure are in nice 
qualitative agreement with what we see on Fig. 6 presenting experimental data on 
e'^e~ annihilation and r decays. 

In summary, our model spectral density (0), with the appropriate values of 
parameters, properly captures all important features which must be inherent to 
spectral densities in real QCD. First, the corresponding polarization operator is 
a sum of an infinite number of simple poles on the unphysical sheet, so that the 
correct analytical properties ensue automatically. The expansion has the right 
structure. The spectral density following from Eqs. ( ^U] ) and ( pUD is dynamically 
smeared by the resonance widths. Purely pictorially it closely resembles what is 
measured experimentally. 



5 Conclusions 

In this paper we address the issue of the preasymptotic component of the spectral 
density not seen in practical OPE. This component oscillates, with the amplitude 
being exponentially damped. The approach to the problem is complementary to 
that of Refs. 0, ^- It is gratifying to observe that the general features of the 
overall picture come out the same as in the previous analyses. Details are different. 
For instance, the instanton-based model discussed in Refs. [|], yields the interval 
between the successive oscillations growing with s, while our present result implies 
equidistant oscillations in the s scale. In this aspect the instanton-based model 
is seemingly less realistic. Moreover, for obvious reasons it does not allow one to 
trace the proper Nc dependence, while our present analysis does. Needless to say 
that it reproduces the desired regularity - the fact that the exponent governing the 
exponential damping of the oscillations is proportional to 

The important lesson we draw is confirmation of a general pattern of the duality- 
violating component in the spectral densities at high energies inferred previously: 
exponential character, modulated by oscillations. Particular details are model- 
dependent. The 't Hooft model is solvable, and all questions can be explicitely 
answered. In four-dimensional QCD it is possible to provide only educated guesses. 
Further efforts are needed to back them up by more solid calculations. 

We believe that the suggested ansatzfoi the spectral density, Eqs. (^) and (pO|), 
gives a very good idea of the duality- violating contributions. It is compatible with 
all general principles of field theory, 1/Nc expansion of QCD, and folklore knowledge 
which is universally believed to be true. 
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Although our discussion was phrased in terms of the spectral densities, its im- 
plications are wider. In particular, it is quite probable that a component of the 
very same structure is present in the so called hard quantities without OPE, e.g. 
thrust. Recently it was realized that the perturbative predictions for such quanti- 
ties must be supplemented by 1/Q corrections (for a review see e.g. [^). If we are 
right, in the intermediate domain of moderate momentum transfers the oscillating 
component might be noticeable too. 

An interesting question which deserves a new dedicated analysis (independently 
of our model of duality-violating contributions) is the behavior of the resonance 
widths as a function of the excitation number at high n in multicolor four- dimensional 
QCD. Another obvious direction of expansion is combining the ip function ansatz 
(^) with information on specific low-dimension condensates which might allow one 
to obtain a fully realistic description of the polarization operator in the entire com- 
plex plane. 
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